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Abstract 

We study the Euclidean two-point function of Fermi fields in the S'C/(2)-Thirring 
model on the whole distance (energy) scale. We perform perturbative and renormalization 
group analyses to obtain the short-distance asymptotics, and numerically evaluate the 
long-distance behavior by using the form factor expansion. Our results illustrate the use 
of bosonization and conformal perturbation theory in the renormalization group analysis 
of a fermionic theory, and numerically confirm the validity of the form factor expansion in 
the case of the £'£7(2)-Thirring model. 
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1. Introduction 

Correlation functions in 2D integrable models have attracted much attention from ex- 
perts in Quantum Field Theory (QFT). The possibility of their exploration on the whole 
length (energy) scale is of great importance. It gives a rare opportunity to probe non- 
perturbatively general principles of QFT. There is also a pragmatic reason for this in- 
terest. The past two decades have witnessed experimental work to identify and study 
quasi one-dimensional systems (for a review, see fll]||). There were collective efforts of 
many physicists to apply integrable QFT to describe such physical systems @. For this 
purpose, a non-perturbative treatment of the correlation functions in integrable models 
seems to be valuable. It is worth noting that in recent years, angle resolved photoemission 
spectroscopy has matured into a powerful experimental method for probing the electronic 
Green's functions in quasi one-dimensional systems [|J. Hence two-point fermion correla- 
tors in integrable theories deserve special consideration. 

In this paper we are studying Schwinger's function (Green's function in the Euclidean 
region) 

(* ff (aO$ ff /(0)) 

in the S'?7(2)-Thirring model, which is described by the Euclidean actionB 

a= I d 2 xl ^^9^^ + ^ OH^*) 2 } . (i.i) 

Here \& a is a doublet of Dirac Fermi fields, and the Pauli matrices r = (r , r 2 , r 3 ) act on the 



"colour" indices a = f, |. The QFT ( |TTT|) possesses a variety of interesting properties 0. 
For instance, it is an asymptotically free theory (for g > 0) with unbroken chiral symmetry, 
and its mass scale M appears through dimensional transmutation. The model belongs to 
the very special class of field theories which admit an infinite number of conservation 
laws preventing particle production in scattering processes. The Hamiltonian of (|1 . 1|) was 
diagonalized by the Bethe Ansatz techniques in 

Also, it is a popular model for the interacting one-dimensional electron g as^ as is 
known PJj JlQjm , ^ , C P- - H) describes the scaling limit of the half-filled Hubbard chain, 

+oo , 111 

H-Hub = - Yl \ ( C L C i+l^ + C Ul,aCj,a ) + U ( <* A C jt T - - ) ( c) A C jA - - J \ . 

j=-oo ^ cr=1,l 



1 The definition of the coupling constant g is not conventional, but convenient for our purposes. 
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where { cj- CT , Cj> >cr i } = $0-0-' t^y . More precisely, if one sends the coupling constant U — > +0, 
the correlation length 

R c = — = e c 
2^ 

diverges and the correlation functions in the Hubbard chain at large lattice separations 
assume certain scaling forms. In particular, if \j— j'\ ^> 1, the equal-time fermion correlator 
can be written as 



0)^)^8^ ^|r~^ F(\f-j\/R c ). (1.2) 



The scaling function F here is directly related to the field-theoretic correlation function: 

<<M*) *.(0) ) = S -^^F(M\x\) . (1.3) 

Z7T \X\ Z 

Our analysis of the fermion correlator is based, on the one hand, on recently proposed 
expressions for the form factors of soliton-creating operators (or topologically charged 
fields) in the sine-Gordon model [|T3[P, and on the other hand, on a conformal pertur- 
bative analysis of two-point correlation functions involving such fields. The form factor 
expressions can be used to obtain the long-distance behavior of these two-point functions, 
whereas Conformal Perturbation Theory (CPT) gives their short-distance expansion [16]. 



The interest in some of these topological fields stems from their role in fermionic theories. 
For instance, it is well-known that the sine-Gordon model is equivalent to the massive 



Thirring model ||17|| . The components of the Thirring fermion field are then associated 
with soliton-creating operators of topological charge ±1 and Lorenz spin ±=, and corre- 
lators of these operators in the sine- Gordon model are related to fermion correlators in 



the massive Thirring model ||18|| . More interestingly, the sine- Gordon theory is closely 



related to a model which is an integrable deformation of (|1 . 1|) ||. This "deformed" (or 



anisotropic) £ , £7(2)-Thirring model exhibits the so-called spin-charge separation, which is 
translated by its representation in terms of two bosonic theories, one for the charge part, 
one for the spin part. The spin part of the fermion field corresponds to soliton-creating 
operators of topological charge ±1 and Lorenz spin ±| in the sine-Gordon model, and its 
charge part is related to similar operators in a free massless bosonic theory. 



2 Without taking normalization into consideration, some of such form factors were considered 



previously in Refs. [14,15| 
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Although form factor expansions and CPT are very effective tools for the study of, 
respectively, the long-distance and the short-distance asymptotics of Schwinger's functions 
|T(||T9| , |2T)|| , one usually gets into trouble when trying to compare both predictions in a region 



where they are expected to be accurate enough. Indeed, in general, one has the freedom 
of choosing the overall multiplicative normalization in the CPT expansion as well as in 
the form factor expansion, and there is no systematic way of relating both normalizations. 
For the case of the soliton-creating operators, the constant relating both normalizations 



was conjectured in |13|]. It allows one to make unambiguous numerical predictions on 
the correlation functions of soliton-creating fields on the whole distance scale using the 
combined CPT and form factor data. We performed this calculation for the case of the 
5't/(2)-Thirring fermion. 

The paper is organized as follows. In Section 2, we recall some standard results con- 
cerning the anisotropic S , £/(2)-Thirring model and its relation to the sine-Gordon theory. 
In Section 3, the short-distance behavior of correlators of the soliton-creating operators 
is examined by means of CPT. Here we also perform a Renormalization Group (RG) re- 
summation of the perturbative expansion in the vicinity of the Kosterlitz-Thouless point 
which corresponds to the SU(2) limit of the fermion theory. In Section 4, the perturbative 
calculation is adapted to the momentum space fermion Schwinger's function; we give the 
two-point function in the 5'C7(2)-Thirring model to third order in the running coupling. 
This particular result was recently obtained by standard perturbation theory in the mod- 



ified Minimal Subtraction (MS) scheme |21| (calculations in |H| concern, in fact, fermion 
correlators in a general non-abelian Thirring model). We then compare the result of Ref. 



21f| with ours and explicitly relate our RG scheme to the MS scheme. In Section 5, the 



long-distance behavior of the fermion correlator in the anisotropic ££7(2)-Thirring model 



is analyzed by means of the form factors given in [|T3[. In Section 6, we examine properties 
of the fermion spectral density in the 5 , C/(2)-Thirring model. The outcome of our calcula- 
tions is discussed in Section 7, where we numerically compare the short-distance behavior 
of the scaling function F ( |1.2| ), (|1.3| ) (from the RG analysis) with its long-distance behavior 
(from its form factor expansion). 

Finally, we note that this work has an essential overlap with Ref. |]22|| , where ( |1 . 1] ) 
was considered as a model of one-dimensional Mott insulators. 
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2. Bosonization of the anisotropic 5£/(2)-Thirring model 

The 5£/(2)-invariant Thirring model admits an integrable generalization such that the 
underlying 577(2) symmetry is explicitly broken down to U(l) <S> Z^: 

Adtm = fd 2 x\ ^ ^7,^ + m + ^ ( JlJl + jpi ) } , (2.1) 
where 

= (2-2) 

are vector currents. The model ( |2.1| ) is renormalizable, and its coupling constants g\\, g± 
should be understood as "running" ones. In particular, in the RG-invariant domain g» > 
\g±\, all RG trajectories originate from the line g± = of UV stable fixed points, and (|2.1j ) 
indeed defines a quantum field theoryi. Hence, in this domain (which is the only one that 
we discuss here), each RG trajectory is uniquely characterized by the limiting value 

p=ilim0||(*) (2.3) 

of the running coupling gn(£) at extremely short distances {£ stands for the length scale), 
i.e. the theory ( |2.1| ) depends only on the dimensionless parameter p, besides the mass 
scale M appearing through dimensional transmutation. 

As is well known (see e.g. 0,§J), the model ( |2.1| ) can be bosonized in terms of the 
sine- Gordon field <p(x), 

A sG = j d 2 x I ^- (d v <p) 2 - 2^ cos(/3vp) | , (2.4) 

with the coupling constant (3 in ( |2.4} ) related to p ( f2.3| ) by 

P 2 = T^- , (2-5) 
1 + p 

and a free massless boson. Then the mass scale M is identified with the mass of the 
sine-Gordon solitons, which is related to the parameter p by [531 



' v^r(i + i; 

M — 2 t 2p 

2r(i) 



2p 
1 + P 

(2.6) 



3 The Hamiltonians corresponding to opposite choices of the sign of g± are unitary equivalent, 
so the sign of this coupling does not affect the physical observables. 
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The precise operator relations between ( |2.1| ) and (|2.3| ) can be found in ]13|. In particular, 
for the two-point fermion correlator, the bosonization implies that 



<tt ff (s)^(0)> = ^ ^ F$(r), (2.7) 
Zn \ x \2 

where we use the notation (n = 1, u = 1/4) for the real function which depends 

only on the distance r = \x\ (and implicitly on the the mass scale M and the parameter 
p), and which, in essence, coincides with the Euclidean correlator of nonlocal topologically 
charged fields in the model ( [2.4| ): 

/ 7 \ tan 

(0-^)0^(0)) = (e w -) #(r), (2.8) 

where z = x 1 + zx 2 , z = x 1 — zx 2 . Again we refer the reader to the paper [1^] for the precise 
definition of the field O™ (a = uj/3). Here we note that it carries an integer topological 
charge n, a scale dimension 

and a Lorentz spin ion. 



3. Short-distance expansion 

3.1. Conformal perturbation theory 

We now turn to the analysis of the short-distance behavior of the correlator ( j2.8| ) . In 
general, one can examine this behavior via the operator product expansion, for instance: 

F™(r) = d(r) + C cos(/M (r) (cos(^) > + .... (3.1) 

The structure functions (Ci(r), C cos ( ( 3 y) )(r), etc.) admit power series expansions in p 2 , 
which can be obtained by using the standard rules of CPT, whereas the vacuum expec- 
tation values of the associated operators are in general non-analytical at p = 0. In the 
perturbative treatment, we regard the sine-Gordon model (|2.4| ) as a Gaussian conformal 
field theory 

•Aoauss = J d 2 x ^ (d„ip) 2 (3.2) 
5 



perturbed by the relevant operator cos(/3(p). Notice that in the limit \x — > 0, the nonlocal 
topologically charged fields O™ can be expressed in terms of the right and left moving parts 
of a free massless field ip = (Pr{z) — <Pl(z) governed by the action 



O:\^ - 6: = exp | i ( a - ± ) <p R (z) - i ( a + ^ ) ^(z) | . (3.3) 
CPT gives the structure function Ci ( |3.1[ ) in the form 

(e- -Y n Ci(r) = (d-^x)Ol uj(3 (0)exp fa f d 2 y cos(^)) \ , (3.4) 

\ Z / \ \ J J I Gauss 

where ( . . . ) Gauss is the expectation value in the Gaussian theory Acauss and the exponen- 
tial is understood as a perturbative series in \i. In the perturbative series, the integrals will 
have 'power law IR divergences which should be thrown away ||16|| . Such a regularization 
prescription is indicated by the prime near the integral symbol. In the absence of loga- 
rithmic divergences, throwing away the divergences is equivalent to treating the integrals 
as analytical continuations in the field dimensions |T^]. Considering only the part of 
perturbative in n, it is a simple matter to obtain 



fW(r) = r~ 2d jl + J n {2uj(3\ -2(3 2 ) ^ 2 r 4 " 4 ^ + o(r 8 " 8 ^, r 2 ) j , 



(3.5) 



where d is given by (|2.9|) and 



J n (a,c)= I d 2 xd 2 y x a+ %x a ~% {l-x)- a -%(l-x)- a+ %x 

y— ty-«+5 (1 - y) a +t (1 - y) a ^ \x - y\ 2c . 

Two comments are in order here. First, the next omitted term in the short distance expan- 
sion ( P^D comes from either the next term in the perturbative series for Ci (0(r 8 " 8 ^ )) or 
from the leading contribution of cos(/3<^) (0(r 2 )) in ( |3.1|) . Therefore, the fj 2 term written 
in (|3.5| ) is a leading correction to the scale invariant part of the correlation function for 
| < (3 2 < 1 only. Second, in writing ( |3.5| ) we specify the overall multiplicative normaliza- 
tion of the nonlocal topologically charged field 0™^ by the condition 

F {n)^^ r -2d ag r ^ _ ( 37 ) 



(3.6) 



The integral (|3.6| ) can be calculated using, for instance, techniques illustrated in ||24j| . 
The result can be expressed in terms of two generalized hypergeometric functions at unity: 

A(q> c ) = 3^2(-c, -c - 1, 1 - q; -c - q, 2; 1) 

(3.8) 

B(q, c) = 3 F 2 (q, q + 1, c + 2; c + q + 2, c + q + 3; 1) . 
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With q = a + n/2 and q = a — n/2, we found: 



J n (a, c) = + + J(3) + j(4) ^ ( 3. 9) 

where 

J (1) =qq T(l - g)r(l - q) T(l + c + q)T(l + c + q)T 2 (-l - c) x 
( cos(7t(q — q)) — cos(nc) cos(n(q + q + c))) A(q, c)A(q, c) , 

J(2) ^\r(l + cni + q W + c + qn 1-c-q) 

T(q)T(-q)T(3 + c + q) w ' ; w ' ' ' 

J(3) ^gr(i + c)r(i + g )r(i + c + g)r(-i- c - g ) 

J(4) _ _ ^(1 + q)TH2 + c)r(-l - c - g-)r(-2 - c - q) 



Notice that for n = 0, the integral ( |3.6|) was calculated previously in the work |25j (see 
also Ref. 



3.2. Renormalization group resummation 

Here we discuss the short-distance expansion of the correlator (|2.8| ) for (3 2 sufficiently 



close to unity For this purpose, it is convenient to use the notation 

e = l-/3 2 <l. (3.10) 
Our previous CPT analysis suggests the following expansion for the structure function Ci: 



Ci(r)=r 2d {l + f; Cfc (^) 2fc }, (3.11) 
^ fe=i J 



where the coefficients are given by certain 4/c-fold Coulomb-type integrals. Evidently, 
this expansion cannot be directly applied in the limit e —>■ 0, where the perturbation 
cos (/?</>) of the Gaussian action (|3.2|) becomes marginal. However, being expressed as a 
function of the scaling distance Mr, the structure function Ci(r) should admit the following 
form: 

d(r)= Cj rm) (Mr) , (3.12) 

where the r-independent renormalization constant Z n , w absorbs all divergences at e = 
and renders the renormalized structure function Cj r finite in this limit. The divergences 



of the renormalization constant Z n>w should be directly related to the singularities of Cj 
at Mr = 0; they point out that the power law asymptotic behavior ( |3.5| ) is modified by 
logarithmic corrections at e = 0. 

In order to explore the short-distance behavior for e <C 1, it is convenient to return 
to the fermion description. Being essentially the corresponding structure function in the 
renormalizable QFT ( j2.1| ), Ci(r) obeys the Callan-Symanzik equation. Therefore it can 
be written in the form: 

C,(r) = r~ 2d exp | - 2 jf ^ (T g - d) J . (3.13) 

Here the function F g is supposed to have a regular power series expansion in terms of the 
running coupling constants g\^± = <7||,_l(?"): 

oo 

r 9 = £ -rtk a\ gf ■ (3.14) 

l,k=0 

Notice that only even powers of the coupling g± appear in this expansion (see footnote 
#3). In writing ( |3.13j ), we use the normalization condition ( |3.7| ), and take into account 
that the U V limiting value of T g coincides with the scale dimension ( |2.9|) , 



\imT g =d. (3.15) 

r— »0 

We have also assumed that there is no resonance mixing of the operator 0™^ with other 
fields, so it is renormalized as a singlet. One can easily check that this is indeed the case 
for the operators with \u\ < \ + 

Condition ( |3.15| ) already encloses an important restriction on the series ( ^.14j ) . Indeed, 



using Eqs. (|2.3|) and ( |2.9|) along with the condition that the line of UV stable fixed points 
corresponds to g± = 0, one obtains 

T 3 = r^u) + T^( g]l )gl + T^( gi] )gj_ + 0(g 6 ± ) , (3.16) 

where 
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The values of the other coefficients 7z j fc>i appearing in G3.14Q essentially depend on the 



choice of a renormalization scheme, i.e on the precise specification of the running coupling 
constants. The latter obey the RG equations 

r d n _ 9± 



dr f\\(9\\,9±) (317) 

dg± 9\\9± 

r = . 

dr f±(g\\,g±) 

Perturbatively, f\\(g\\,g±) and f±(g\\,g±) admit loop expansions as power series in g» and 
g±. In this work, we will use the scheme introduced by ALB. Zamolodchickov |2^,^7|]. He 
showed that under a suitable diffeomorphism in g» and g±, the functions /n and f± can 
be chosen to be equal to each other, and furthermore, to be equal to 

/||=/± = l + f- (3-18) 

With this choice for the /3-function, the RG equations ( |3.17| ) can be integrated. To do this, 
we note that this system of differential equations has a first integral, the numerical value 
of which is determined through the condition 



g\-gl = {2p) 2 . (3.19) 
Using (|3.19|) , (|3.10|) and (|2.5|), equations (|3.17|) are solved as 



9 1| = 2 P > 9± = P t^— , 3.20 

11 1 — q 1 — q 

where 

1 — n\ ~ 2e 

-y-j =(rA) 4e . (3.21) 

The normalization scale A is another integration constant of the system ( |3.17| ). It is of the 
order of the physical mass scale and supposed to have a regular loop expansion, 

A = M exp (r + Tip + r 2 p 2 + . . .) . (3.22) 

It should be noted that the even coefficients To, T2, ■ ■ ■ are essentially ambiguous and can be 
chosen at will. A variation of these coefficients corresponds to a smooth redefinition of the 
coupling constants which does not affect the /3-function. By contrast, the odd constants 
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r 2 k+i are unambiguous and precisely specified once the form of the RG equations is fixed. 
It is possible to show []23| , p7| that the odd constants vanish in the Zamolodchikov's scheme: 



r 2k+l = (k = 0, 1 . . .) . 
Once the coefficients r 2 k in ( P-22| ) are chosen, the running coupling constants are 



completely specified, and all coefficients in the power series expansion ( 3.14 ) are determined 
unambiguously. They can be explicitly calculated by comparing the CPT result (|3.5| ) with 
the form ( gig) . From ( gl3| ), 

V 9 = -\ ^rlog (Ci) 



and, as it follows from the general CPT expansion (|3. 11 ) and the definition ( |3.21|) of q 



the function T g can be expanded in powers of q. Explicitly, using the CPT result (|3.5|), 
T g = d-2e (^) 4f V Jn(2w(l-e),2e-2) q + 0(q 2 ) . (3.23) 

Moreover, the coefficients in this expansion are power series in p. For example, using Eqs. 
( FDD and ( |3.22|) , it is easy to show that 



T (^A ) 2£ = exp { " 2toP+ ( 2ro " I) p2 ~ ( 2r2 + 2r ° " I C(3) " I) p3 + 0(/) } • (3 ' 24) 



Here and after, we set for convenience 



e ro = y | e 7E+ "° , (3.25) 
where = 0.5772 ... is the Euler constant. The integral J n (2u(l — e), 2e — 2) appearing 



in (|3.23|) can also be expanded in powers of p, using e = p/(l + p). In Appendix A, we 
quote the first few terms in the expansion of J n (a, c) ( |3.6| ) around c = —2, which are 
obtained through the use of (|3.9| ). From this expansion, it is easy to obtain the expansion 
of J n {2u{l — e), 2e — 2) in powers p. Then, one can compare the CPT expansion of F g 
in q and p (|3.23|) with the corresponding expansion (|3.14| ) coming from the RG analysis 
(where of course one should expand grn and g\ in q and p from ( |3.20| )). This determines 
the coefficients 7^1 for I = 0, 1,2. If we want an expression valid to order gr 4 , we need 
one more coefficient: 70,2- In principle, it can be obtained from the expansion in p of the 
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coefficients in the series ( 3.11 ). In Section 5, we describe a way to find 70,2 without the 
cumbersome calculation beyond the lowest CPT order. 



In order to simplify the form of the structure function (|3.13| ), it is convenient, instead of 
using the coefficients 7/^, to parametrize the first few terms of the power series expansions 
rC 1 ' 2 )^) (OD as: 



r ( %ii) 



1 



n u\ 

i + f i32-y +Wl ^ + ^ 2 

V2 



(3.26) 



r (2) (^,i) = -f + o( 9ll ). 



The explicit values of the coefficients iti, U2, v± and v-i in ( |3.26j ) are given in Appendix B. 

Let us substitute (|3TTBD and flOp into Eq. ( gTg ). The RG flow equations QCT] ) 
allow one to evaluate the integral and to write the structure function in the form ( |3.12j ) 
with 

-4^ 2 -n 2 (l+p 2 )/4 I 2 ^ 2 -n 2 (l-p 2 )/16 

X 

■" ltf " U2S " ^l+^(^i+^||) + 0(/) 



Ci ren) =(Mr) 



(3.27) 



and 



M 2d [2 p+1 ^e T0p+T2p3+ 



o \ n 2 /2-2d . 

' ' e 2pu 1 +(2p) 3 u 2 + . 



(3.28) 



Notice that the transformation 



e w +w 1 (2p) 2 +w 2 (2p) 4 - + ... z 



(3.29) 



where the series contains only even powers of p with arbitrary coefficients Wk, accompanied 
by the transformation 



(ren) 



{ A 2 \ { A A \A 1 ( 

e -wo-w 1 (g ll -g ± )-w 2 {g {l -g ± ) +... 



(ren) 



does not affect the structure function Ci ( |3.12j ) due to relation ( |3.19p . 

Our prime interest in this work is the correlation function ( |2.7| ). For n = 1 and u = j, 
the relations obtained above lead to the following perturbative expansion for the two-point 
fermion correlator in the anisotropic S , £7(2)-Thirring model: 



(^(x)# a (0)> 



2tv 



lpf__ 

2+4 



\X 



(9l) - exp \ - n - T JLg\ 



16 



32 



x 



In 2 3 



-]> i - yr n -J ,1 - - (To - g r - + 0(, 4 ) 



(3.30) 
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where 



„3 



n \ f 3p lE 3 4 , 



^ = (4p)-TO (MJ- I e xp^-^V + 0(p 



In Eq. ( |3.30|) , we use the notation tq defined by ( |3.25|) . 



We now set p = and gii = g± = g in ( |3.3U| ) to obtain the perturbative expansion of 
the scaling function F ( |1.3|) for the 5't/(2)-Thirring model, 



F ^ = exp | _ ± g + JL ( f0 _ 1 ) 9 , _ JL ( f . _ i.) g 3 + 0{gi) } . (3. 31 ) 

Here the running coupling constant g solves the equation 

-^-i + Iln(^) = ln(J ( |e^+^Mr) , (3.32) 

which is the limit p = of Eqs. (ggTj) and (|372lD . 

Let us stress here that, if the perturbation series could be summed, then the function 
F should not depend on the auxiliary parameter Tq\ 

^ = °- 

This is, however, not true if we truncate the series (|3.31[) at some order iV (for instance, if 
one leaves only the terms explicitly written in ( |3.31|) ). In this case, 

9 ^(pert) 



8t 



F% ert) = 0(g N+1 ) 



where the truncated series is denoted by F^ ert \ In fitting numerical data with (|3.31| ), we 
may treat Tq as an optimization parameter, allowing us to minimize or at least develop a 
feeling for the effects of the remainder of the series. Similar ideas have been discussed for 
QCD in Ref. |8 . 



It may be worth mentioning that Eq. (|3.27| ), along with explicit values of the coefficients 
quoted in Appendix B, allows one to immediately determine the short-distance expansion 
of some other conventional correlators in the (anisotropic) 5 , [/(2)-Thirring model. For 
example, since the sine-Gordon field ip ( p.4j) itself can be defined by the relation 
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n = 
a = 



and the spin current (|2.2|) is bosonized as 



4 = ( 3 - 33 ) 



we can use ( |3.27| ) to obtain the short-distance expansion of the current-current correlator. 
For the £'£7(2)-Thirring model ( |1 . 1| ) one has, 

( J fa) 4f (0) ) = { ( <W " ^w) h + <W h) , (3.34) 



where 



^ = l-| + (^ + i)f-^o + 4 + (f + ^ + ^-^-^C(3))/ + 0( 9 =) 
/ 2 = ^ {l - 2r„ S + f„(3f„ + - (4fS + I fS + | - H _ * ((3) ) g> + Ofe*)} , 



and the running coupling constant g is the same as in (3.32). 



4. Perturbative expansion of the momentum-space correlation function 

4-1. Large-momentum asymptotics 

Perturbative calculations of fermion Green's functions in renormalizable 2D models 
with four-fermion interaction are widely covered in the literature. The results in this 
domain are usually expressed in momentum space. Hence it seems appropriate at this 
point to adapt the calculation of the previous section to the momentum-space fermion 
correlator, giving a large-momentum expansion. 

The RG analysis performed in the previous section can be applied in essentially the 
same way to the Fourier transform of the fermion correlator ( [2.71) : 

j d 2 xe~ ipx (v^(x)^(O)) = -zcW ^ F(p 2 ) . (4.1) 

Here and after we use the notation p 2 = p^p^. From the result of CPT, (|3.5| ), one imme- 
diately obtains the large momentum expansion of this Fourier transform: 

F = Q(d 9 ) (p 2 ) d *- l > {l + Q( q*~J €) ~2(3 2 ) ^V)" 2e + 0({p 2 )-*\p- 2 ) }, 

(4.2) 
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where 

Q(a) = 2 1 ~ 2a r( 2" fl ) 



r(i + o) 

and 

, 1 ^ 
= - + 



2 4(1 + p) 

is the scale dimension of the fermion field. The factor Q(d^ — 2e) /Q(d^) is essentially the 
only source of differences between the RG treatments in coordinate space and in momen- 
tum space. The RG analysis in momentum space goes as in the previous section. The 
perturbative part in p of F obeys the Callan-Symanzik equation, so it can be written as 

pipert) = Q{d ^ {p 2 )d ,-\ exp |_jf^ (f g - d V ) } , (4.3) 

where the function F g admits a power series expansion in terms of the momentum-space 
running coupling constants g»± = g\\ : ±(p 2 ) depending on the Lorentz invariant p 2 : 

oo 

r, = E %* 9\\9f ■ ( 4 - 4 ) 

l,k=0 

Notice that, with some abuse of notations, we use here the same symbols g»± for the 
momentum-space running couplings as we used for the coordinate-space running couplings. 
In order to fix the coefficients in ( |4.4| ), we have to choose a renormalization scheme. Sub- 
stituting r by 1/ \pp* in ( |3.21| ) defines Zamolodchikov's scheme in momentum space. It is 
a simple matter to repeat the steps of the previous section in order to determine the first 
few coefficients 7^1 in (|4.4|) . Just compare the logarithmic derivatives of the expressions 
( [Op and (fO|); the only difference is that the factor Q(d^ — 2e)/Q(d^) in (|4.2|) will have to 
be expanded in p, giving non-trivial contributions. As for the coefficients 7^2, one would 
in principle need the next order in CPT. However, again as in the previous section, it is 
possible to determine 70,2 without this calculation, as described in the next section. From 
these coefficients, and from the form of the RG flow equation, one can evaluate the inte- 
gral in ( |4.3| ) and obtain the asymptotic behavior of the two-point function in the Euclidean 
region at p 2 — > +00. We quote here the result in the case of the £ , £7(2)-Thirring model, 

= exp { - 1 9 + 1 ( f0 - i) t - A [rl - 1) i + Otf) } . (4.5) 

Here 

-g- 1 + - ln(g) = ln(^Me fo /\/^) , (4-6) 
14 



and To is an arbitrary parameter which can be chosen at will. Notice the strong similarity 
between QDJ and fl£3ip . 

We also quote here the corresponding function T g ( |4.4| ) in the case grn = 

1.3 2 3_ a , 3 ( „ _9 , _ 3 



4-2. Comparison with the four-loop conventional perturbation calculations 



In [211, the anomalous dimension for the fermion field in the MS scheme was found to 



fourth order for a general non-abelian Thirring model (see also [£9] and references therein 
for a discussion of various aspects of dimensional regularization in the non-abelian Thirring 
model and for results to lower order). In contrast, we have calculated, in coordinate space, 
the two-point functions of more general operators, including the fermion fields, in the 
particular case of the S'C/(2)-Thirring model (and an anisotropic deformation of it), and 
we have sketched the equivalent calculation in momentum space for the fermion fields. We 



would now like to compare Eq.( |4.5| ) with the SU(2) case of the Ali-Gracey result |2T 



In order to perform the comparison, we need to find the relation between our running 
coupling constant g and theirs, which will be denoted qag = — ^ ^ an d then find the 
relation between our function F g (O) and their anomalous dimension, which we will 
denote 7a- 

The coupling A corresponds to the MS scheme; the associated /^-function was found 
in |2]| to fourth order: 

V ^2 = & = + £ ~ A ' + °( A ") • ( 48 ) 

ap z 7T 2n z 1287I- 3 

By comparison, in the scheme that we use, the (3- function ( |3.17| ), ( |3.18| ) is 

2 " 2 ^ = "» = -T^72 = -« 2 + T-T + ^- ^ 

The difference in the factor multiplying the square of the coupling in these two expressions 
results only from a different normalization of the coupling in the action (see Eq.( |l.l|) ). The 
relation between the couplings g and A that corresponds to these different /3-functions is 

A rgZ+(T* + ^ + ^-)g3 + 0(g*). (4.10) 



7T ° J V 2 128 



4 Notice that in pl[| , the coupling constant gAG is assumed to be negative, so A > 0, which 
agrees with the sign of our coupling constant g. 
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Here r is some numerical factor which cannot be determined by comparing the /3-functions: 
its variation modifies the choice of the normalization scale and doesn't affect the /3-beta 



functions. The normalization scale for the MS scheme is defined by imposing the following 
condition on the subleading asymptotics of the solution of the RG flow equation (|4.8|) : 

A _ 1 1 \n\n(y^/A m ) | Q fin 2 ln( y^/A^ 

n ln(^/A m ) 2 ln 2 (VF/A^) 



(4.11) 



(This implies that the term 0(1/ In (y^/Aj^g-)) does not appear in the expansion of A.) 
From flU), ( CT ) and (fCTl|), we find that 



2tvM e 



TQ—T 



(4.12) 



In ]2T| , the perturbative part of the function F ( |4.1| ) was calculated up to the overall 
multiplicative normalization to third order in A. The result can be written in the following 
form 



1 i 






oc - — exp < 
h\ 1 







where the function h\ and the anomalous dimension 7> were given in |2T| to fourth order 
in A for the Thirring model with a general non-abelian symmetry. In the particular case 
of the 5'L r (2)-symmetry, they specialize to 



hy 



15 ^2 
1 + „ A 2 



11 



1287T 2 



512tT 



A^ + 



3(80C(3)-511) 



and 



7a 



32768 7T 4 
3 



A 4 + 0(A 5 ), 



A + 0{\ l 



(4.13) 



(4.14) 



16tt 2 64tt 3 1024tt 4 

Comparing ( f4.3|) in the case p = with the above expressions, one has the following 
relation: 

9 2 2 2 dX BK 1 



Using Eqs. ( f4.10|) -( fl.l4|) , one can check that our result ([O]) agrees with (|4.15|) , provided 
that 

(4.16) 



T = T . 

Notice that the relation between the normalization scale A^jg and M, 



A 



MS 



2tvM 



(4.17) 



which is a consequence of ( f4.12| ) and ( |4.16|) , was previously found in Ref. (12 
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5. Long-distance behavior 

Here we concentrate on the long-distance behavior of Schwinger's function ( |2.8| ) for 
n = 1 and u = 1/4. Let us recall that for | < (3 2 < 1, there are only solitons and 
antisolitons in the spectrum of the sine-Gordon model. We will denote them by A_ and 
A + respectively. The conservation of the topological charge, 

p r a « 

— / dxd x v? , 

^ J-oo 

implies that the non-vanishing form factors of the operator Opj 4 are of the form 

( vac | O+} 4 (0) | A-{0 X ) ■ ■ ■ A.(6 N+1 ) A + (9[) ■ ■ ■ A+(9' N ) ) , (5.1) 

where $i and 9j denote rapidities of solitons and antisolitons respectively. Up to an overall 
normalization, all these form factors can be written down in closed form, as certain iV-fold 



integrals [5D|,ST],H|. The spectral decomposition for the correlation function ( |2.7| ) then 
gives 



d6 l de 2 de^ 

■x 



. /-l-oo m -i pc 

?» = L S e ~ (e, l<"-l o ?»l^-C))l 2 + l/_ 0o (ar) . 
e - Mr EL, J2 K« ac I °t/M I A n (e 1 )A a2 {e 2 )A,M))\ 2 + --- , 

cri+(T2+cr 3 = — l 

(5.2) 

where the dots stand for the five-particle and higher contributions, which are of the order 
of e _5Mr . The long-distance asymptotic behavior of the correlation function is dominated 
by the contribution of the one-particle states, 

( vac | O+} 4 (0) | A_{6) ) = ^ (e+if) , 

and has an especially simple form, 

( p — Mr 

FWW = Z 1 (/3/4)| 7 = + 0( e —)}. (5.3) 



Here we use the notation Z n (a) (a = u/3) from work [13] for the field-strength renormal- 
ization which controls the long-distance asymptotics of the correlation function (|2.8| ) . Let 
us stress here that the overall multiplicative normalization of the field was already 
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fixed by the condition (3.7|), hence the constant Zi(/3/4) is totally unambiguous. In [13 
the following explicit formula for Z n (w/3) was proposed: 



z„M) 



Co 



si n 

Lo \ 4 



2C\) \16pJ 



0FMr(l + fa 



2r(i + 4 



2d 



exp 



cosh(4wi) e -( 1+ ^ nt - 1 



+ 



n 



t I 2 sinh(t) sinh ( (1 + p)t) cosh(tp) 2 sinh(t) 



2d e 



-2t 



(5.4) 



In this formula, d is the scale dimension (|2.9|) and the constants Ci, C2 read 



Ci = 



12 



ir(|) 



exp 



dt sinh 2 (f)e- t 
* 2 cosh 2 (tp) sinh(t) 



Co 



r 4 (|) 

47T 3 



exp 



2 - 



dt sinh 2 (f)e" 



t cosh(tp) sinh(t) 



where Aq = 1.282427. . . is the Glaisher constant. 

We do not write down explicitly the general formula for the three-particle contribution 



in ( |5.3|) because it is a rather mechanical substitution of relations presented in [13]. (For 
f3 2 = 1 the corresponding formulas can be found in Appendix C.) Here we make the 
following observation concerning the /3 2 — > 1 limit. The examination of ( |5.2j ) based on 
explicit formulas for the form factors shows that the function [Zi(/5/4)] 1 admits 
an asymptotic power series expansion in terms of the variable p 2 . In other words, all 
divergences at p 2 — > of Fy^ considered as a function of the variables p 2 and Mr, are 
absorbed by the normalization constant Zi(/3/4). Using Eq. (|5.4|) , one can check that 
the constant Zi(/3/4) admits exactly the same type of singular behavior at p 2 = as the 
constant Z n ^ ( |3.28| ) for n = 1, u> = 1/4, and also that 



= 2-i ^e"* Aq exp ( Wl p 2 + 0(p 3 )) 

Z l,l/4 



(5.5) 



The explicit form of the coefficient w\ is not essential here. What is important is that 
the linear term in p does not appear in the expansion ( |5.5|) . This observation can be 
immediately generalized and checked for any n and ui. Furthermore, we expect that 



log 



Z n (u/3) 



J2^k p 2k +o( P c 



(5.6) 



where Z n (w/3) is the normalization constant ( |5.4j ). In other words, by means of the trans- 
formation ( |3.29| ) with properly chosen coefficients w^, the constant Z n ^ in (|3.12| ) can be 
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set to be equal (in a sense of formal power series) to Z n (ujf3). At the moment, we do not 
have a rigorous proof of (|5.6|) . But it leads to some interesting prediction to be checked. 
As was already mentioned, the calculations performed in the leading CPT order determine 
only the combination v<i — 3u2/2, but do not fix the individual values of the coefficients 
U2 and t>2 in the series ( |3 . 2 7| ) . Accepting (|5]6|), one can immediately find the values of 
the coefficients U2 (see Appendix B). In the case n = 1, u> = |, it allows one to extend 
the perturbative expansion ( |3.31| ), as well as the equivalent expansion Ql.5| ), to order g 3 . 
As was discussed in Section 4, Eq. ( f4.5| ) is in a complete agreement with the result of 
four-loop perturbative calculations from [^TJ . This in fact shows that the p 3 -term really is 
absent in the series (|5.5|). 



6. Spectral density 

The spectral density is an important quantity related to the two-point function and its 
analytical structure in momentum space. It is often what is measured in actual condensed 
matter experiments [f|,|22| , and it allows one to completely reconstruct the two-point func- 
tion. In this section, we discuss the properties of the spectral density in the 5'C7(2)-Thirring 
model. 

The spectral decomposition of the fermion Green's function yields the following form 
for the function F ( f4.1|) : 

AF(s) 



F{p ) = 1 ~ / 2 ds ^T7 • (fU) 

Jm 2 P + s 

The notation AF for the spectral density reminds us that F, considered as a function of 
one complex variable p 2 , has a branch cut in the Minkowski region p 2 < starting at 
p 2 = — M 2 , and that the spectral density can be recovered from the discontinuity along 
this cut: 

AF(s) = ^- {F{e™s)-F(e-™s)) . (6.2) 

The easiest way to obtain the large s asymptotics of the spectral density is to use the 
expansion ( f4.5|) along with knowledge of the analytical properties of the coupling constant 
g as a function of the complex variable p 2 . Notice that g can be expressed in terms 
of the principal branch of the product log (or Lambert) function, which gives the solution 
for W in W e w = z (see e.g. @): 
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The principal branch of the ^-function analytically maps the complex z-plane minus the 
branch cut z £] — oo, — e _1 ] to the part of the complex VF-plane enclosing the real axis and 
delimited by the curve Jfe W = -3m W cot (Sm W) for — tx < Sm W < tc. The analyticity 
implies that the power series 

n=0 

converges for real positive s > e~ l and \(j)\ < it and coincides with W(e l< ^s). Similar 
considerations are, of course, valid for the coupling constant g In particular, for 

sufficiently large s, 



00 1 A 



. . dp z 

This then gives us, with ( }4.5| ) and the RG flow equation (|4.9|) , the asymptotic expansion 
of the spectral density for large s. It can be written in the following form: 



A«.>~£{l-§-^ + (*}^ 



(6.4) 



Here the function F^ pert ^ is given by (|4.5[) and (7 is defined by the equation (fOJ). 



Now let us consider the threshold behavior of the spectral density. According to 
the analyses of the previous section, the long-distance asymptotic behavior of the scaling 
function F ( |i~3| ) is described by the expansion 

F = F (1) +F (3) +0(e" 5Mr ) , (6.5) 



where 



with the constant 



F (l) = C e -Mr 



C = 2 - 6 e -i A% = 0.921862 . 



The function in O gives the three-particle contribution to the correlation function. 
Using the definitions (|1.3|), (fO) and the above relation, one can obtain: 
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Here the dots stand for contributions of the massive multiparticle intermediate states. The 
last relation implies that the spectral density (|6.2|) can be written as 

7T y's/M 2 - 1 

where 

1 for s > 



^ I for s < 

and AF( 3 ) is some function which contributes to the spectral density only above the 
threshold s = 9M 2 . 



7. Numerics 

In Table 1 we present results of numerical evaluation of the function F ( |1.3|) as a 
function of the scaling distance Mr (r = |x|). To estimate the short-distance behavior, we 
use the perturbative expansion ( |3.31| ). As was already mentioned, the parameter fo allows 



one to have control over the accuracy of the truncated series, so we calculate ( |3.31| ) for 
two different values of fo : —0.25 and +0.25. To determine the long-distance behavior of 
the function F, we use the formula (|6.5|) , where the three-particle contribution F^ s ' was 
obtained by means of Eq. ( |5.2| ) along with formulas for the three-particle form factors 
quoted in [[13| (see Appendix C). It is interesting to see that the sum of the one- and three- 



particle contributions to F is very near to unity at r = (to within 1%), which indicates 
that this three-particle computation of the correlation function is in fact accurate to about 
1% for all distance scales (more accurate, of course, for larger r). Also, note that the 
crossover between the long- and short-distance asymptotics appears to be at the scaling 
distances Mr ~ 0.001 — 0.01, where both asymptotics coincide to within about 0.1%. 
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Mr 


pW 


i?(3) 




F (pert) ( fQ = _ _25) 


F (pert) ( fQ = o 25) 





.921862 


.068 


.990 


1.00000 


1.00000 


.00001 


.921853 


.0553 


.9771 


.980129 


.980130 


.00005 


.921816 


.0522 


.9740 


.976311 


.976314 


.0001 


.921770 


.0504 


.9722 


.974192 


.974196 


.0002 


.921678 


.0483 


.9700 


.971674 


.971678 


.001 


.920941 


.0415 


.9624 


.963508 


.963520 


.002 


.920020 


.0375 


.9575 


.958435 


.958454 


.01 


.912689 


.0252 


.9379 


.939386 


.939460 


.025 


.899101 


.0168 


.9159 


.919294 


.919494 


.05 


.876902 


.0106 


.8875 


.894050 


.894547 


.075 


.855251 


.00738 


.86263 


.871796 


.872717 


.1 


.834135 


.00541 


.83955 


.850520 


.852013 


.15 


.793454 


.00317 


.79662 


.808380 


.811548 


.2 


.754757 


.00200 


.75676 


.765139 


.770842 


.25 


.717947 


.00131 


.71926 


.719980 


.729252 


.3 


.682932 


.000889 


.683822 


.672640 


.686654 


.35 


.649625 


.000617 


.650243 


.623153 


.643171 


.4 


.617942 


.000436 


.618379 


.571774 


.599063 


.45 


.587805 


.000313 


.588118 


.518942 


.554677 


.5 


.559137 


.000227 


.559365 


.465257 


.510405 



Table 1. The scaling function F ( |1.2| ), ( |1.3| ). The first columns give the 
results of the long-distance expansion which includes contributions of the one-, 
three- and one+three-particle states. The data in the last two columns correspond 
to the perturbative expansion ( |3.31|) for the two different values of the auxiliary 
parameter tq. 



Appendix A. 

In this appendix, we give the first few terms in the expansion of J n (a, c) ( |3.6|) around 
c = —2. The coefficients in this expansion involve standard functions of a, which could 
then easily be used to obtain an expansion of J n (-pp-, jqr~) i n powers of p, as is needed in 
( p.23| ). To simplify the result, we will use the parameter 

6 = c + 2. 

We find the following expansions in 6 of the functions A(g, 6 — 2), B(q, 6 — 2) ( |3.8| ) involved 
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in (HOD: 



T{2-b)T(l + b-q) [ V 6 
Hence, 

t , r r\ \ 71-2 4a 2 - n 2 + 2nb f . . _ 

J n (a, 6-2) =- _ ^ 2 — exp j - G n {a) b+ 

( G 'n(a) » oGl(a) + G n (o) 10 \ 3 
I 12 2 n2-4«2 + 3 C ^ j r + ( ^ j 

where 

G n (a) = ip(a + n/2) + ^(-a + n/2) + 2^ E , 
and G' n (a) = £G n (a), G«(a) = £ G n (a) . 



Appendix B. 

In this appendix, we write down explicit expressions for the coefficients u±, u 2 , v± and 
V2 taking part in the expansion ( }5.16| ), ( pj.26| ) of the the function T g . 

On the one hand, from the assumption ( |5.6| ), the coefficients of odd powers of p in the 
exponential factor of Z n>UJ ( |3.28j ) are completely fixed by the conjectured constant Z n (uf3) 
(|5.4|) . This fixes u\,u 2 uniquely, giving 

Ul ^4) (Tn(2w) _|) + ^, 

n(n + 4) llu 2 1 T2 / o n 2 \ 

+ 777^0 + "77 I" - 777 



768 48 ' 24 ' w ' 2 V" 16 ) ' 

(B.l) 

where 

T n (o) = V(o + n/2) + ^(-o + n/2) + 2 lE + 2f , 

d d 2 
T M = -T- T n (a) , T n (a) = T n (o) . 

On the other hand, the expansion in powers of p of ( |3.23| ) uniquely determines the 
coefficients ui, v\ and v 2 — § «2 in the first equation of ( |3.26| ). The coefficient u\ thus 

23 



obtained is in agreement with ( |B.1|) , verifying the assumption ( |5.6|) to first order. The 
coefficients v\ and (using the expressions from (|B.1|)) i> 2 are 



16/2 4 V 16 



-i))u<(2a;) 



1 / o n 2 \ . 1 / 9 n(n — 4)\ ~. , 1 / 9 n(n — 2) 1\ m . 

1/.2 ^ (o i A sfo\ o\ n(n-8) U! Vl 3u 2 f 
-8\ U "16 J (2T2 - 14C(3) - 3) -^56— + Y + Y + — • 

Appendix C. 

In this appendix, we give the formula for the three-particle contribution ( |6.5|) to 
the fermion two-point function in the 5'C/(2)-Thirring model that we used for our numerical 



calculations. We first specialize the expression written in |L3| to the case of three-particle 

0/4 



form factors of the field 0\ /A for (3 2 = 1: 



vac I Ol /4 (0) I A_{9x) . . . A+(6 k ) . . . A_(6 3 ) ) in = ^ x 

3 /c 3 

f[e^ J[G{e m -e j )\ j ^e-i[[^- 7 ) n ^(7-^p)+ (C.l) 

7i=l m<j < JC+ p=1 p=k+l 

jc - % p=i p=k * 



p=L p 

Here the functions G and W are 



G(9 ) = , 2 ljfi£tij si „h( 9 /2) exp f f * ^r^f' ) (C.2) 



and 



y/nA% k ' j \J t sinh(2t) cosh(t) 



W3 _ ifl \p/ 1 1 i9_ 



WW = 2 V4 ^ 2( ^ 4 ^ . (C.3) 

The contour C+ starts from —00 on the real axis of the complex 7-plane, goes above the 
poles located at 7 = 9 P + in/2, p = 1, . . . , k and below those located at 7 = 6 p — in/2, 



24 



p = k + 1, . . . , 3, always staying in the strip — tt/2 — < Sm7 < tv/2 + 0, and finally 
extends to +oo on the real axis. Similarly, the contour C_ goes above the poles located 
at 7 = 9 P + m/2, p = 1, . . . , k — 1 and below those at 7 = 9 P — iir/2, p = k, . . . , 3. Notice 
that the integrals in (|C.1|) can be expressed in terms of the generalized hypergeometric 



function 3F2 at unity. 

Using the expressions ( |C.1[ ) and performing one of the rapidity integrals in ( |5.2| ) , one 
can obtain the following form for the function in 



2g — r°° -Mr^/3+2 cosh x+2 cosh y+2 cosh(x-y) 

_p(3) _ G / rfxdy X 

3yrr 6 (i) ' (3 + 2coshx + 2coshy + 2cosh(x-?/))i 

1. 

(2\R 1 (x,y)\ 2 + \R 2 (x,y)\ 2 ) \G(x)G(y)G(x - y)\ 2 ( - l±^± ^ ' 

y e -+- e y -+- 1 

The functions R± and R2 here are 

R2(x,y) = e"^ +i ? R^-x,y- x) - e~^~^ R\(-y,x - y) 

and 

cosh I cosh I tj ( 1 1 ix 1 iy zx iy 



Ri(x,y) 2sinhxsinhy 2' 2 2tt' 2 2tt' 2tt' 2tt 

« cosh cosh § /l 1 i(y — x) 1 ix „ i(y — x) „ ice 

+ e~ 2 ? ? — 7T( i — I • 1 1^ i 2 H 

2sinh(y - x) sinhx V2'2 2tt '2 2tt' 2tt ' 2n 

-I cosh ^ cosh f /l 1 _ z( x - y) 1 _ _ z(x - y) zy 

6 2sinh(x-y)sinhy V2'2 2tt '2 2tt' 2tt ' 2tt 

where C/(a, 6, c; d, e) is related to the generalized hypergeometric function 3 F 2 by 

r(a)r(6)r( c ) 

Z7(a, fr, c; d, e) = ph\f73 3-F2K 0, c; d, e; 1) . 
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